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Abstract. 

Using the international ground-based network of two-frequency receivers of the 
GPS navigation system provides a means of carrying out a global, continuous and 
fully-computerized monitoring of phase fluctuations of signals from satellite-borne 
radio engineering systems caused by the Earth's inhomogeneous and nonstationary 
ionosphere. We found that during major geomagnetic storms, the errors of 
determination of the range, frequency Doppler shift and angles of arrival of 
transionospheric radio signals exceeds the one for magnetically quiet days by 
one order of magnitude as a minimum. This can be the cause of performance 
degradation of current satellite radio engineering navigation, communication and 
radar systems as well as of superlong-baseline radio interferometry systems. 
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1. Introduction 

Radio engineering satellite systems (RESS), with their ground-based and 
space-borne support facilities, are finding ever-widening application in various 
spheres of human activity. They are able to provide global coverage, accuracy, 
continuity, high reliability and meet a number of other requirements imposed 
when tackling a broad spectrum of engineering problems. However, the use 
of RESS also implies new (and, in some cases, more stringent) requirements 
dictated by the need to ensure safety and economical efficiency of the operation 
of ground-based and airborne facilities, as well as to solve special problems 
(observation, aerophotography, searching and rescue of distressed transport vehicles 
and people). This applies equally for performance of Global Navigation Satellite 
Systems (GNSS) as well as for very-long-baseline radio interferometers (VLBI) 
[Thompson, R. et al, 1986]. 

Such requirements cannot be met unless the influence of destabilizing factors 
in the radio wave propagation channel is taken into account. As electromagnetic 
waves propagate through the ionosphere, they experience quite varied disturbances 
[Davies, 1969; Kolosov et al, 1969; Yakovlev, 1985; Goodman and Aarons, 1990; 
Yakubov, 1997]. 

The key characteristic of the ionosphere that determines the variation of radio 
wave parameters is the integral (total) electron content (TEC) I(t) or its derivatives 
(with respect to time and space) l' t , l' x and i' along the propagation path. 

TEC variations may be arbitrarily classified as regular and irregular. Regular 
changes (seasonal, diurnal) - for the magnetically quiet mid-latitude ionosphere at 
least - are described by models providing relative accuracy of TEC prediction in 
the range 50. ..80%. Irregular changes (variations) are associated with ionospheric 
irregularities of a different nature, the spectrum of which has a power law character 
[Gajlit, 1983; Gershman et al, 1984; Yakubov, 1997]. 

TEC variations introduce proportionate changes of the signal phase 
ip(t,x,y) = kil(t,x,y), which gives rise to measuring errors of the range 
aD = k 2 dl, the frequency Doppler shift of the signal af = k 3 l' t , and the angles 
of arrival of the radio wave aa x = k^l' x and aa y = k/J.' y , because the last 
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four quantities are proportional to the time and space derivatives of the phase. 
Furthermore, the maximum value of the measuring error of angular deviations can 



The coefficients k\ — k 4 are inversely proportional to the signal frequency f c or 
to its square [Davies, 1969; Kravtsov et al, 1983; Goodman and Aarons, 1990]. A 
calculation uses a Cartesian topocentric coordinate system with the axis x pointing 
eastward E, and the axis y pointing northward N. 

Investigations of phase fluctuations of transionospheric signals have been and 
are carried out using radio beacons on satellites with circular and geostationary 
orbits [Komrakov and Skrebkova, 1980; Livingston et al, 1981; Gajlit, 1983]. The 
trouble with these measurements is that temporal and spatial resolution is low, and 
continuity and global coverage of observations are unavailable. 

The use of the international ground-based network of two-frequency receivers 
of the navigation GPS system that at the beginning of 2002 consisted of no less 
than 1000 sites and is posting its data on the Internet, opens up new avenues for a 
global, continuous, fully computerized monitoring of phase fluctuations of signals 
and associated errors of RESS performance. 

Some research results on the prediction and estimation of radio signal 
fluctuations and errors of RESS performance caused by them were reported in 
earlier work [Afraimovich and Karachenschev, 2002]. 

Results of experimental investigations of phase fluctuations of RESS signals 
that are discussed in this paper were derived from analyzing the data for a set 
consisting of 100 to 300 GPS stations covering the time interval 1998-2001, for 
different conditions of geomagnetic activity (the D st -index from -13 to -377 nT). 
The number of passes of GPS satellites with a duration of no less than 2.5 hours, 
the data from which were used in the analysis, totaled no less than 300000, or 
several orders of magnitude in excess of the currently available statistic of such 
measurements. 

Below we give an outline of the techniques used in this study and illustrate 
their application in the analysis of ionospheric effects. The overall sample statistic 
of errors oD, af and era is presented for different geomagnetic conditions. To ease 



be deduced using the relation era = /c 4 J(I' X ) 2 + (-C) 2 - 




5 



comparison with other research results reported in [Gajlit, 1983; Kravtsov et ai, 
1983], the errors <jD, of and aa are calculated for the working frequency of 300 
MHz. 

2. General information about the database used 

This study relies on the data from the global network of receiving GPS stations 
available on the Internet (Fig. 1). As is evident from Fig. 1, the receiving sites are 
relatively dense on the territory of North America and Europe, and less as dense in 
Asia. Fewer stations are located on the Pacific and Atlantic. 

Such coverage of the terrestrial surface by GPS receivers makes it possible, 
already at the present time, to address the problem of a global investigation 
of ionospheric disturbances and their consequences with a very large spatial 
accumulation. 

Thus, in the Western hemisphere the corresponding number of stations is 
as large as 500, and the number of line-of-sight (LOS) to the satellite is at least 
2000. ..3000. This provides a number of statistically independent series at least two 
orders of magnitude higher than would be realized by recording VHF radio signals 
from first-generation geostationary satellites or low-orbit navigation satellites - 
TRANSIT [Gershman et al, 1984]. 

Our analysis used the North-American sector with a large number of GPS 
station - range of longitudes -120° ... -60° E and of latitudes 20° ... 70° N. Table 1 
presents information about number of the days, the number n of GPS arrays 
composed of three stations, the data from which are used, and extreme values of 
D stmin . Total number of spectra E also is resulted which were used for definition of 
average sizes of inclinations. 

For a variety of reasons, for the various events to be analyzed, we 
selected somewhat differing sets of receiving stations, yet the geometry of 
the experiment was virtually identical for all events. The coordinates of the 
receiving stations used in the experiment were available at the electronic address: 
ftp: / /lox. ucsd.edu/pub.processing/gamit/setup/coords.txt . 
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3. Analysis of the measuring errors of the range, Doppler 
frequency and angles of arrival of the radio wave caused 
by changes in the regular ionosphere 

Recently a number of authors [Wilson et al, 1995; Mannucci et al, 1998; 
Schaer et al., 1998; and others] have developed a new technology for constructing 
Global Ionospheric Maps (GIM) of TEC using IONEX data from the international 
IGS-GPS network. The GIM technology and its uses have been reported in a large 
number of publications [Wilson et al., 1995; Mannucci et al., 1998]. 

The standard IONEX format is described in detail in [Schaer et al, 1998]. 
Therefore, we will not give a detailed description of the GIM technology for 
reasons of space but limit ourselves only to the information required for the 
presentation of our method. Two-hour TEC maps are easily accessible to 
any user, which are calculated by several research groups in the USA and 
Europe and are availabe on the Internet in the standard IONEX format 
([ftp: / / cddisa.gsfc.nasa.gov /pub/ gps / products / ionexj ) . It is also possible to obtain 
15-min maps if necessary. 

Fig. 2 is a schematic representation of a single elementary GIM cell. The cell 
nodes are designated as a, b, c, d. The cell size (5° in longitude and 2.5° in latitude) 
is determined by the IONEX file standard. For simplifying the transformations to 
an approximation sufficient for our problem for latitudes not exceeding 60°, the cell 
can be represented as a rectangle with the sides d e and d n . It is easy to overcome 
this limitation by complicating to a certain extent the transformations allowing for 
the sphericity; however, we do not present them in this report. 

The linear size of the rectangular cell in latitude is independent of the latitude 
and is d e = 279 km; the linear size in longitude depends on the latitude, and for 
40°iV it is d n = 436 km. 

For each time t, for the nodes a, b, c, d from the IONEX file the values of 
vertical TEC are known - I a , I^, I c , Id. 

The determination of the range D by the phase method is based on measuring 
the phase difference p between the received signal and the reference signal formed 
in the receiver. Such a measurement can be made at the intermediate or carrier 
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frequency of the signal. In this case: 

where c is the propagation velocity of radio waves in a free space. 
Generally the quantity ip for the transionospheric propagation may be regarded 
as the sum of two components [Afraimovich et al, 1998]: 

<p = y? s +A</9 (2) 

where ip s is the main component associated with a change of the distance 
between the signal source and the receiver. 

Analysis of the measuring errors of the range, Doppler frequency and angles of 
arrival of the radio wave caused by changes in the regular ionosphere investigations 
of global phase variations of radio signals and their influence on the operation of 
RESS ought to take into account the proportionate relationship between phase 
(phase derivative) changes of the transionospheric signal and corresponding TEC 
variations [Kravtsov et al, 1983; Goodman and Aarons, 1990]: 

A(p = 8.44 x 1(T 7 x I -^f^ + v9 (3) 

Jc 

where f c is the radio wave frequency (Hz); I a (b,c,d) is the TEC measured at the 
points a, b, c, and d (10 16 el / m 2); and ip is the initial phase [Spoelstra and Kelder, 
1984]. 

By way of example we now analyze the errors of measurement of the range 
aD, the frequency Doppler shift af and the angle of arrival of the radio wave a a 
using the IONEX data and the method that was developed at the ISTP SB RAS 
[Afraimovich and Karachenschev, 2002]. 

The change of the range that is introduced by the ionosphere (ionospheric 
error) is proportional to Aip: 

aD = cx (4) 

Upon substituting (3) into (4), we can obtain the expression for determining 
the measuring error of the range aD introduced by the ionosphere: 
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c x 8.44 x 1(T 7 xdl 
aD = = 4.48 x dl (5) 

2 n Jc 

As is seen from (5), the error aD is directly proportional to the TEC variation 
dl and inversely proportional to the carrier frequency squared. 

Using the values of the spatial derivatives of TEC l' x and I' y and of the 
derivative of TEC with respect to time l' t makes it possible to uniquely obtain - for 
each instant of time - the values of errors of determination of the angle of arrival 
era and the frequency Doppler shift af by formulas [Kravtsov et al, 1983; Goodman 
and Aarons, 1990]: 

1.39 x 10 2 



= — J 2 x V^) 2 + W ( 6 ) 

1.34 x 1(T 7 _, 
aa = x I t (7) 

Jc 

In the simplest case, the values of the derivatives for the selected cell of the 
map can be obtained using TEC increments for the four cell nodes and for two 
times t\ and t 2 — t\ + d t : 

dl = (I a2 - Ll + h2 ~ hi + Ic2 - hi + Id2 ~ /dl)/4 

Al' t = (I a2 - I al + I b2 - I bl + I c2 - I cl + I d2 - I dl )/4d t 

y°) 

AI X = (J cl - I bl + I dl - I al + I c2 - I b2 + I d2 - I a2 )/Ad e 
Al' y = (I al - I bl + I dl - I cl + I a2 - I b2 + I d2 - I c2 )/4d n 

Where necessary, the spatial derivatives can be estimated by taking into 
account the TEC values in adjacent nodes of the map, and the time derivative 
(with a time resolution of IONEX maps no worse than 15 min) can be inferred by 
averaging increments for several successive time counts. 

The procedures (5), (6) and (7) are performed for all cells of the selected 
spatial range and for the selected time interval of the day. One variant of data 
representation implies a full exploitation of the IONEX format with the difference 
that, rather than the values of TEC and the error of TEC determination [Schaer et 
al, 1998], are entered into the corresponding cells of the map values of error aD, 
af and aa. 
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By way of example, it is appropriate to give the results derived from analyzing 
the regular part of the spatial-temporal TEC variations for a relatively magnetically 
quiet day of July 29, 1999 (with the largest deviation of the D st - index of -40 nT) 
and for the magnetically disturbed day of April 6, 2000 (with the largest deviation 
of the IVindex of -293 nT). 

Fig. 3 a, b, c portrays the maps of the errors aD, af and aa obtained on the 
basis of files in the IONEX format for the magnetically quiet day of July 29, 1999 
in the geographic coordinate system in the range of longitudes -120° ... -60° E and 
latitudes 20° ... 70° N . Fig. 3 d, e, /, respectively, characterizes the values of aD, 
af and aa for the magnetically disturbed day of April 6, 2000. The figure also 
shows the time interval 19-21 UT, for which the analysis was carried out. Contours 
show the values of errors of phase measurements in units, respectively, of aD - 
"m" (meters), af - "Hz" (Hertz), and aa - "arcmin" (minutes of arc). The vertical 
calibrated scale shows the maximum and minimum values of the corresponding 
errors. 

The above maps are a pictorial rendition of the behavior dynamics of the 
errors aD, af and aa in the spatial and temporal ranges selected. Noteworthy is a 
considerable difference of the maps for the magnetically quiet and disturbed days. 
As is evident from Fig. 3, gradients of spatial distribution of the errors during 
disturbances increase more than an order of magnitude when compared to the quiet 
period, which would lead to a degradation of RESS performance. 

4. Analysis of the irregular errors aD, af and aa 

Our analysis of the irregular component of the errors was based on using raw 
data in the form of series of TEC values for selected receiving sites as well as values 
of elevations Q s (t) and azimuths a s (t) to visible satellites that were calculated using 
a specially developed program, CONVTEC, to convert RINEX-files (standard files 
for the GPS system) available from the Internet [Afraimovich et al, 1998]. 

Fig. 4 presents the geometry of transionospheric radio sounding. The axes 
z, y, x are pointing, respectively, to the zenith, the north N, and to the east E; 
P is the point of intersection of the LOS to the satellite (the line connecting the 
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satellite to the radio signal receiver) with the ionospheric F 2 region peak; S is the 
subionospheric point (projection of the point P onto the terrestrial surface). 

The GPS technology provides the means of estimating TEC variations on the 
basis of phase measurements of TEC / in each of the spaced two-frequency GPS 
receivers using the formula [Hofmann-Wellenhof et al, 1992]: 



J = ah one J 2 r 2 [(Ll\l ~ L 2 \ 2 ) + Const + TlL], (9) 



J_ IV? 

40.308 fl - /| 

where L\\\ and L 2 \ 2 are phase path increments of the radio signal, caused 
by the phase delay in the ionosphere (m); L\ and L 2 are the number of full phase 
rotations, and Ai and A2 are the wavelengths (m) for the frequencies f\ and f 2 , 
respectively; const is some unknown initial phase path (m); and nL is the error 
in determination of the phase path (m). 

Series of the values of elevations Q s (t) and azimuths a s (t) of the beam to the 
satellite were used to determine the coordinates of subionospheric points, and to 
convert the "oblique" TEC I (t) to the corresponding value of the "vertical" TEC 
I it) by employing the technique reported by [Klobuchar, 1986]: 



I = In X COS 



urcsin ( — 1 cos©. 



(10) 



where R z is the Earth's radius, and h max —300 km is the height of the F 2 -layer 
maximum. 

All results in this study were obtained for elevations Q s (t) larger than 30°. 

To analyze the errors aD, erf and era that are cased by the irregular component 
of TEC variation, we make use of the relations for the respective regular errors (5), 
(6) and (7). The difference in the analysis of irregular errors implies a different 
(compared with the IONEX technique) method of determining the values of TEC 
and its derivatives [Afraimovich and Karachenschev, 2002]. 

The phase differences A<p X;y along the axes x and y are proportional to the 
values of the horizontal components of TEC gradient G E — l' x and G E = l' y - 

To calculate the components of the TEC gradient V x and I' y are used linear 
transformations of the differences of the values of the filtered TEC (Ib — I a) and 
(Ib — Ic) at the receiving points A, B and C: 
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j, = Va{Ib ~ Ic) ~ Vc{Ib - I a) , j, = x c (Ib - I a) - x a (Ib - Ic) ^ 
xaVc ~ x c Va ' V %aVc ~ x c Va 

where xa, Da, xq, Uc are the coordinates of the sites A and C in the topocentric 
coordinate system. When deriving (11) we took into account that xb = Ub = 0, 
since site B is the center of topocentric coordinate system. 

The time derivative of TEC V t is determined by differentiating I(t) at the 
point B. 

The procedures of (5), (6) and (7) are performed for all groups of three GPS 
stations of the selected spatial range and for satellites visible from these stations, 
as well as for the selected time interval of the day. 

Fig. 5 presents the results derived from analyzing the irregular component of 
the errors aD, erf and era in the form of fluctuation spectra of the range, Doppler 
frequency and angles of arrival of radio waves. 

With the purpose of improving the statistical reliability of the data, we used the 
spatial averaging technique for spectra within the framework of a novel technology 
[Afraimovich et al, 2001]. The method implies using an appropriate processing of 
TEC variations that are determined from the GPS data, simultaneously for the 
entire set of GPS satellites (as many as 5-10 satellites) "visible" during a given 
time interval, at all stations of the global GPS network used in the analysis. 

Individual spectra of the errors aD, erf and era were obtained by processing 
continuous series of I(t) measurements of a duration no shorter than 2.5 hours. 
To eliminate errors caused by the regular ionosphere, as well as trends introduced 
the motion of satellites, we used the procedure of removing the linear trend by 
preliminarily smoothing the initial series with the selected time window of a 
duration of about 60 min. 

Fig. 5 shows the overall character of the TEC variations dl(t) that were filtered 
from TEC series obtained from measurements of the phase difference between two 
coherently coupled signals from the GPS system [Hofmann- Wellenhof 1992] for the 
magnetically quiet day of July 29, 1999 (panel a, at the left) and a major magnetic 
storm of July 15, 2000 (panel e, at the right). Furthermore, the panels show the 
station names and locations, as well as GPS satellite numbers (PRN). 
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As is evident from the figure, the intensity dl(t) during the disturbance even 
at such low latitudes is increased an order of magnitude as a minimum. This is 
reflected on logarithmic amplitude spectra IgS(F) of TEC fluctuations and their 
derivatives (left-hand scale in the figures) and of the fluctuations of aD, af and 
era, converted to the working frequency of 300 MHz (right-hand scale) which are 
represented on a logarithmic scale (panels b, c, d, f , g, h). 

The logarithmic amplitude spectrum IgS(F) obtained by using a standard 
FFT procedure. Incoherent summation of the partial mplitude spectra lgS(F)i of 
different LOS was performed by the formula: 



lg(S(F)) = lg 



E S(F)i 

i=l 

n 



(12) 



where i is the number of LOS; i — 1, 2, ... n. 

As a consequence of the statistical independence of partial spectra, the 
signal/noise ratio, when the average spectrum is calculated, increases due to 
incoherent accumulation at least by a factor of ^fri, where n is the number of LOS. 

Fluctuation spectra from the magnetically quiet day of July 29, 1999 are 
shown by the thin line in panels /, g, h (Fig. 5) for comparison with the spectra 
from the disturbed day. The range of fluctuation periods is shown in bold type 
along the abscissa axis in panels d and h. Panels b and / show also the number n of 
GPS arrays composed of three stations, the data from which are used to estimate 
the spatial derivatives of TEC [Afraimovich et al, 2001]. 

Spectra of phase fluctuations have a power law character with the values of 
the slopes v, shown in panels b, c, d, f, g, h. The slope of spectrum is determined 
by the slope of the fitted straight line (thin black line in panel b of Fig. 5). These 
results are in reasonably good agreement with data reported in [Komrakov and 
Skrebkova, 1980; Livingston et al, 1981; Gajlit et al, 1983; Kravtsov et al, 1983; 
Gershman et al, 1984; Yakubov, 1997]. 
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5. Distribution of slopes and scales of spectra of the errors 

In the preceding section we obtained the fluctuation spectra of the errors aD, 
af and aa corresponding to different ionospheric conditions. It was shown that 
the spectra have a power law character, with definite indices of the slope v (Fig. 5). 
However, the individual spectra that were obtained do not give a full insight into 
the global picture of fluctuations of the errors. To obtain a generalized estimate 
we carried out an analysis of the spectra of errors obtained by processing the data 
covering more than 30 days (over 600 spectra). 

In our experiment the range of fluctuation periods of spectra varies from 2 min 
to 2 hours. Experimental spectra of the errors aD, af and aa have a complicated 
form (Fig. 5). There are spectral density maxima and minima. Generally, however, 
the spectral density S(F) decreases with the increasing fluctuation frequency. For 
each individual spectrum, calculated on a logarithmic scale, we determined the 
fitted straight line. 

The result of a generalized estimation is represented by the plots in Fig. 6. 
The plots show the distributions of the slopes P aD (u), P a f(y), P aa (v) of power law 
spectra of the errors aD, af and aa, respectively. The interval Au of the slopes 
of the spectra aD is A aD = —1.45. ..2. 45, with the mean value of <A aD >= —1.95. 
In view of the relationship between the spectrum aD ~ <p(t) and the spectra 
af and aa that are proportional to the derivative of (p(t,x,y) (af ~ ip '(£) and 
aa ~ <p'(x, y)), for the slopes of the spectra af and aa these values are, respectively, 
A af = -0.40.. .1.60, <A af >= -0.99 and A aa = -0.40.. .2.00, <A aa >= -1.8. 
Our results are in good agreement with findings reported by a number of authors 
[Gajlit et al, 1983]. 

The fitted straight line (the slant thin black line in panel b, Fig. 5) for the 
logarithmic scale of errors in range may be described by the expression (13): 

Y = axX + b (13) 

here a characterizes the slope of the straight line, and b is the scale coefficient 
characterizing the rise of the straight line with respect to the abscissa axis, that is, 
the value of F when X = 0. 
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The value of lgS aD (F) can be determined from a single (previously determined) 
set of values of lg(F), v aD and b a r> on the basis of the expression: 

lgS° D (F) = v aD xlg(F) + b aD (14) 

Using (14) and the resulting value of b a r> and v a n we can determine the value 
of any spectral component of lgS aD (F) and determine its contribution to the total 
error. Proceeding in a similar way, we can describe the spectral characteristics of 
af and era. 

Panels a, b, c in Fig. 7 presents the distributions of the values of the scales of 
power law spectra of the errors b aD , b a f and b aa obtained for the entire set of data 
(over 600 spectra). 

Let us estimate the amplitude IgS^ of the 32-minute harmonic F of the 
spectrum of the errors aD (Fig. 5) by making use of the values of <v aD >= —1.96 
(Fig. 6) and < b aD >= —6.59 (Fig. 7). Upon substituting the corresponding values 
into (14), we obtain lgSQ D (F ) = -0.153, which corresponds to aD = 0.70 [m]. 

6. Conclusion 

Main results of this study are as follows: 

1. In this paper we have suggested a new technique for estimating the errors of 
RESS performance, based on using phase measurements from two-frequency 
receivers of the satellite navigation GPS system. 

2. Fluctuation spectra of the errors are generally quite well approximated by 
the power function. Mean values of the slopes of spectra, obtained from the 
experiment, are in good agreement with data reported by other authors. 

3. The errors aD, af and era are conveniently estimates with a sufficient 
accuracy for practical purposes by using mean statistical data for indices of 
slopes of spectra v and scale coefficients b. 

4. Our results can be equally assigned to the estimation of the Russian GNSS - 
GLONASS performance quality. 
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Figure 1. The map showing the locations of the receiving stations of the global GPS 
system. The number of receiving sites totals no less than 1000 as of the beginning of 
2002. 

Figure 2. Geometry of an elementary GIM cell. Nodes of the cell are denoted a, b, c, 
and d. The cell size is determined by the values of d n = 279 km and of d e dependent on 
the cell's location latitude. 

Figure 3. Maps of global distribution of measurement errors of the range aD, frequency 
Doppler shift erf, and of angles of arrival of radio waves era (left - for a magnetically 
quiet day, right - for a magnetically disturbed day). The errors can be evaluated from 
the graduated scale, shown in the figure at the right. 

Figure 4. Schematic representation of the transionospheric sounding geometry. The 
axes z, y and x are directed, respectively, zenithward, northward (N) and eastward (E). 
P - point of intersection of Line-of-Sight (LOS) to the satellite with a horizontal plane 
at the height of the maximum of the ionospheric F 2 - region h max ; S - sub ionospheric 
point; and O s , a s - azimuth and elevation of the direction r along LOS to the satellite. 
Figure 5. TEC variations dl(t) for the magnetically quiet day of July 29, 1999 (left, 
panel a) and a major magnetic storm of July 15, 2000 (right, panel e). These panels also 
show the names and coordinates of the COSA and BLYT stations and GPS satellite 
numbers (PRN). Panels b,c, d, f, g, h on a logarithmic scale present the amplitude 
spectra S(F) of TEC fluctuations and derivatives (left scale in the panels), and the 
fluctuation spectra of aD, erf and era converted to the working frequency of 300 MHz 
(right panel). For comparison, spectra from the magnetically quiet day of July 29, 1999 
are shown in panels f,g,hbj the thin line. The range of fluctuation periods is shown 
in bold type along the abscissa axis of panels d and h. Panels b, e show the number n 
of GPS arrays. The power law character of the spectra is determined by the values of 
slopes v. The thin black slant line in panel b is a fitted straight line of the power law 
spectrum. 

Figure 6. Probability distributions Paoiy)-, Pafiy) and P aa (y) of the slopes u a n, v a j 
and v aa of the power law spectra of the errors aD, af and aa. Mean values of the slopes 
v, and also total numbers of the values of the slopes of the corresponding spectra use in 
the analysis are shown in all panels. 
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Figure 7. Probability distributions P b (aD), P&(cr/) and Pb(ua) of the scales Ig b aD , 
Ig b a f and Ig b aa of the power law spectra of the errors. Mean values of the slopes 
< Ig b a o > {^g baf >,< Ig b aa >) (on a logarithmic scale), standard deviations, and also 
total numbers of the values of the scales of the corresponding spectra used in the analysis 
are shown in all panels. 
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